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A VECTOR-DYADIC DEVELOPMENT OF THE EQUATIONS 
OF MOTION FOR N-COUPLED RIGID BODIES 
AND POINT MASSES 


Harold P. Frisch 
Goddard Space Flight Center 


INTRODUCTION 

The design and analysis of most spacecraft attitude-control systems is usually based upon 
the assumption that the spacecraft itself can be modeled adequately as one body or as a 
collection of contiguous rigid bodies. 

When, in the analyst’s judgment, two or more bodies are required for representative 
modeling, the derivation of the system’s equations of motion becomes laborious and 
subject to numerous errors in SIGN and judgment. The general-purpose digital program 
N-BOD has been developed to relieve the analyst of this time-consuming and error-prone 
task. 

N-BOD assumes that the spacecraft can be modeled as a topological tree of rigid bodies, 
momentum wheels, and point masses. It also assumes that: 

• All momentum wheels are symmetric and embedded within rigid bodies and that 
point masses exist only at limb ends. 

• Unless otherwise directed by the user, all gyroscopic interaction torques are 
significant. 

• All contiguous rigid bodies are connected by either zero-, one-, two-, or 
three-axis gimbals and point masses have either zero, one, two, or three degrees 
of relative translational freedom. 

N-BOD provides the user with two modes of output: 

• It may be used to derive and to output in vector-dyadic form on a line printer the 
complete nonlinear equations of motion of the system described to it. 

• It may be used to solve numerically the equations of motion and to output any 
set of system-state variables. 

Several techniques are given in the literature for deriving the equations of motion for an 
arbitrary number of interconnected rigid bodies: namely, References 1 through 6. 
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In each of the formalisms cited, the N-coupled rigid bodies must form a topological tree; 
that is, no closed paths are permitted in the system. 

Two basic approaches are developed in the literature which show how formalistic methods 
can be evolved which will ultimately define the torque-free dynamics of an arbitrary 
N-coupled-body system. 

Hooker and Margulies (Reference 2), and Roberson and Wittenburg (Reference 4) choose 
to write the equation of motion of each individual body of the connected system, taking 
into account in a very formalistic manner the interaction forces between adjacent bodies. 
Velman (Reference 1 ) uses what he refers to as a nested-body approach; that is, rather 
than having N-vector equations, one for each of the N bodies, the bodies ar? grouped into 
N nests and the equation of motion of each nest is then defined. 

While the first method is well documented in the open literature, the nested-body method 
of Velman is outlined only in Reference 1. The equations of the nested bodies can, 
however, be derived from Velman’s discussion and basic principles defined in most texts 
on rigid-body motion; for example, see Reference 7. 

In Reference 3, Hooker comes to the realization that by making use of the nested-body 
concept, the method used by both himself in his first paper and by Velman in Reference 1 
for handling constraints can be vastly simplified. Rather than switch from the discrete- 
body to the nested-body approach as Hooker does in Reference 3. one may attain con- 
siderable simplification simply by starting from the nested-body concept. 

The purpose of this report is to combine the best features of Hooker’s discrete-body 
approach and Velman’s nested-body approach into a computationally efficient formalism. 
Material in this report will concentrate exclusively upon the theoretical development of 
the equations programed in N-BOD. 

BASIC SYSTEM 

The basic system to be studied can be defined as an arbitrary numoer of rigid bodies, 
momentum wheels, and point masses coupled together in such a manner as to form a 
topological tree (no closed paths). 

Contiguous rigid bodies are assumed to be connected together by either a zero-, one-, two-, 
or three-axis gimbal; point masses are assumed to have either zero, one, two, or three 
degrees of relative translational freedom; momentum wheels are assumed to be symmetric 
about their spin axes and to be imbedded within rigid bodies. Rigid bodies which contain 
momentum wheels are referred to as gyrostats. 

Let 

N = total number of rigid bodies, gyrostats, and point masses 
M = total number of imbedded momentum wheels. 
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To define the system mathematically, each body, momentum wheel, point mass, and hinge 
point must be assigned a unique label. The set of consecutive integers will be used. Let 
body 1 be the principal body of the system. All other bodies are given distinct integer 
labels ranging from 2 to N inclusive. For simplicity of computation, it is convenient for 
the labeling of the bodies to be such that, along any topological path from body 1 to its 
end, the body labels are of increasing numerical magnitude. 

As a direct consequence of this labeling convention, the connection topology of the N-body 
system can be uniquely defined by the N X 1 connection matrix J (X), where 

J(l) = 0, and 

J(X) * label of the body contiguous to and inboard of body X, X = 2, 3, . . ., N. 

All hinge points must be assigned unique labels. It is efficient to define them such that: 

Hinge point X - 1 is the point of connection between the contiguous bodies 
J (X) and X. 

All momentum wheels are given distinct integer labels from 1 to M inclusive. These may 
be randomly assigned. 

Let 

MO(m)= body label of the gyrostat in which momentum wheel m is imbedded. 

In order to distinguish between rigid bodies and point masses in the system, two sets of 
body labels are defined: 

S R = set of all body labels associated with rigid bodies, and 

S L = set of all body labels associated with point masses. 

The union of S R and S L is the set of all body labels S; that is, 

S = S R US L . (1) 

To define the equations of motion of the entire N-body system, the system is broken 
up into N distinct nests of bodies. The equation of motion for each nest is then defined 
and solved simultaneously with the equations for all other nests. 

The nests are each given distinct integer labels from 0 to N - 1 , inclusive, such that: 

S k _ ( = set of all body labels associated with those bodies which, relative to body 
1 , are outboard of hinge point k - 1 of body k. The elements of 
S k _, define the bodies which make up the nest k - 1 . 

To define the motion of body X relative to a reference frame fixed at the origin of the 
nest S k l , the bodies lying along the topological path from hinge point k - 1 to the 
center-of-ma's of body X must be known. This information is given by the connection 



matrix J(\); however, it is notationally efficient to define the sets of body labels 
S kixi such that 


S k „ M = set of all body labels associated with those bodies lying on the 
topological path from hinge point k - 1 to the center-of-mass of 
body X. 

As an aid to the visualization of the contents of each set of body labels defined, a 
particular example is given which is general enough to bring out the salient features of 
the notation defined. 

Figure 1 defines a particular 10-body system with each body and hinge pouv given a 
distinct integer label. Note in particular, that the body labels are of increasing numerical 
value along any topological path beginning at hinge point 0; furthermore, momentum 
wheels can be numbered randomly. 



Figure 1. Labeling Scheme for 10-body Example 
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From the figure it can be seen that the contents of each set defined are: 


S R = 

{1,2, 3, 4, 5, 

6, 9,10} 

s L = 

|7,8| 


S 

{ 1,2, 3, 4, 5, 

6,7,8,9,10} 

So = 

{1,2,3, 4, 5, 

6,7,8, 9,10} 

s, - 

{2, 3, 4, 5, 6, 

7,8} 

S, - 

{3,8f 


S 3 * 

{4,5, 6,7} 


S 4 = 

{5,7} 


S 5 = 

{6} 


S 6 • 

{7} 


S, = 

{8} 


S. - 

{9} 


S, = 

{ io} 


MO(l) * 

3 


MO(2) = 

1 


MO(3) = 

5 



5 





The elements of the set S kl X( are given in the following example. 


\x- 1 
k-Y\ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

1 

1.2 

1,2,3 

1,2,4 

1,2, 4, 5 

1,2, 4, 6 

1,2, 4, 5, 7 

1,2,3, 8 

1,9 

1,10 

1 


2 

2,3 

2,4 

2,4,5 

2,4,6 

2, 4,5, 7 

2,3,8 



2 



3 





3,8 



3 




4 

4,5 

4,6 

4,5,7 
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5,7 




5 
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7 
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8 



8 
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All undefined sets in the above table are to be taken as empty sets. Physically, 

S k _, x , = { [implies that body X is not contained in the nest k - 1. 

Each rigid body of the system is assumed to have nonzero mass and rotational inertia. 

Each point mass is assumed to have nonzero mass and zero rotational inertia. 

Let 

m x = total mass of body X plus that of all momentum wheels imbedded in it, 

4> x * »otal inertia tensor of body X plus that of all despun momentum wheels imbedded 
in it, about the composite body center-of-mass, and 

I w = inertia tensor of momentum wheel m about its center-of-mass. 

m 

To define the linear or angular momentum of each body or nest of bodies in the system, 
position vectors locating hinge points and centers of mass must be given. 

Let 

]?, = position vector from origin of inertial reference to hinge point 0 of body 1 , 

= position vector from hinge point J (X) - 1 to the hinge point X - 1 of body X, and 
ii? x = position vector from hinge point X - 1 to the center-of-mass of body X. 




The exact position and labeling of each of these vectors is shown in figure 2 for the 
standard 10-body example. 



Figure 2. Positions and Labeling of Hinge-point and Center-of-mass Location Vectors 

The vectors and dyads appearing in the equations of motion are most naturally defined 
relative to certain body fixed-reference frames. 

Let 

Reference frame 0 = inertially fixed frame of reference having its origin at ;iie 

inertial origin, 

Reference frame X = body X fixed frame of reference having its origin at the hinge 

point X - 1 . If body X is a point mass, the reference frame X 
is taken to be fixed relative to reference frame J(X) and have 
its origin at hinge point X - 1 . 

To avoid explicit definition of coordinate systems and their associated kinematics, all 
equations are written in terms of vectors and dyads. Inevitably, this leads to the necessity 
of differentiating, relative to an inertial reference frame, vectors which are more easily 
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expressed in some other moving frame. Extensive use is made of the identities from 
vector differential calculus that 


+ w, X 


( 2 ) 


where 

V * vector to arbitrary point P from origin of reference frame i, 

Vj = velocity of P with respect to the inertial frame, 

^ = velocity of P with respect to moving reference frame i, 

(3, * angular velocity r»f moving reference frame i with respect to the inertial frame, 

pj * inertial velocity of the origin of reference frame i, 

and 

© ©© 

+ w, X + w, X (w, X v) + 2«, X v + v , (3) 


acceleration of P with respect to the inertial frame, 

inertial acceleration of the origin of reference frame i, 

linear acceleration of P due to angular acceleration of reference 
frame i, 

centripetal acceleration of P due to rotation of reference frame i, 

Coriolis acceleration of P, and 

apparent acceleration of P relative to reference frame i. 


EQUATIONS OF MOTION (VECTOR-DYADIC FORM) 

The entire N-body system is free to translate and rotate relative to a fixed inertial origin. 
All contiguous rigid bodies are free to rotate relative to each other, point masses are free 
to translate relative to their respective contiguous bodies, and all momentum wheels are 
free to rotate about their axis of symmetry relative to the gyrostat in which they are 
imbedded. 


where 

1 - 
Pi 
X? 


c5, X (t3, X t) = 

2(3, x? 

CD© 
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G, . * linear momentum of body X relative to the inertial origin, 

■ t A 

L, x = angular momentum of body X relative to the inertial origin, 

1^ = angular momentum of momentum wheel m relative to the inertial origin, 

F x = resultant of external forces acting on body X, 

- resultant of external torques acting on body X, and 

- resultant of external torques acting on momentum wheel m. 

w m 

Then from Newton’s three fundamental laws of motion, the equations of motion for each 
body X and each momentum wheel imbedded therein are given by 

( 4 ) 

Vv (5) 


H_ * X , 


where 


m: MO(m) = X * all m such that MO(m) * X . 

Similarly, the equations of motion for the nest of bodies k - 1 are 


Lc u -L? x , 

^ L U ' ^ K’ 


H m *^ w » m: M(Xm)e s k-| , 


where 


sum over all bodies X contained in the nest k - 1 , and 
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m:MO(m) eS k ., * all in such that M(Xm) is contained in the nest k - 1 . 

The equations of motion for the entire N-body system are obtained by solving 
simultaneously 2N + M vector equations. These equations may be either in the discrete-body 
form of Equations (4), (5), and (6) or in the nested-body form of Equations (7), (8), and 
(9). This analysis follows the nested-body form. 

Between each pair of connected contiguous bodies, forces and torques of constraint exist 
which limit the number of degrees of relative freedom from six to, at most, three. 

Let 

F' , * resultant force of constraint actLig on body k through hinge point k - 1 , 

, = resultant torque of constraint acting on body k through hinge point k - 1 , and 

= resultant torque of constraint acting on wheel m preventing motion about 
m any axis normal to its spin axis. 

For most problems of practical interest, mechanisms such as springs, dampers, motors, 
and so forth exist between the contiguous bodies. 

Let 

F” j = resultant force acting on body k due to all mechanisms existing between 
bodies J(k) and k at hinge point k - 1, 

0^, = resultant torque acting on body k due to all mechanisms existing betv i 

bodies J(k) and k at hinge point k - 1 , ai:d 

CL m = resultant torque acting on momentum wheel m due to all mech 
existing between it and body M(Xm). 

It should be noted that at every hinge point of the system (with the exception of hinge 
point 0), the forces and torques defined above exist in equal and opposite pairs. Thus, for 
any nest of bodies containing the hinge point k - 1 , these forces and torques are internal to 
the nest and have a resultant effect of zero. At the hinge point of the nest, however, they 
are external. 

External forces may be present which act on the N-body system. These may be locally 
applied, distributed throughout, or applied only over selected portions of the system. (It 
is assumed that external forces are not applied directly to momentum wheels.) 

Let 

= resultant external force applied to point i of body X, and 

R X1 = vector from center-of-mass of body X to the point i at which the force 
is applied. 
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For Newton’s laws of motion to be applied, motion must be defined relative to the 
inertial origin. 

Let 

7|, x * position vector from inertial origin to centerof-mass of body X, and 

7 k l K ■ position vector from hinge point k-1 to centcr-of*mass of body X. 


From the definitions given in *' e previous section, it follows that 



and 






us 




+ O, 


k-l.X-l 
i # k 


00 ) 

(ID 


From these definitions, the resultant of the external forces acting upon the nest k-1 is 
given by 


2 F A - F' , + F» , + 22 F<*>, > 

*Sk., A k *‘ k " 


( 12 ) 


while the resultant of the external torques acting upon the nest k-1 , measured relative to 
the inertial origin, is 



and the resultant torque acting upon momentum wheel m is 




(14) 


The vector equations which completely define the motion of the nest k-1 are obtained by 
direct substitution of Equations (12;, (13), and (14) into Equations (7), (8), and (9). 




♦ FI 


k-l 



(15) 
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Et lk 

x«Q 


♦E, * <?» - ?,) x (Fi- * F«.) 


XfS 


ff u x F« ♦ ♦»>), 


( 16 ) 


kl 


and 


H m = ^ m +CL m 

m 


m: MO(m)e S k j > (17) 

where, for the sake of notation compression, the following definitions have been made: 

p(®) = / 4 F (e) , 
and 

^ 

(19) 


(18) 


f(«) = 


St 


X F<*>. 


ieX 


The inertial momentum for each body X and wheel m may be defined in terms of their 
respective mass and inertia properties and the system-state variables. 


Let 


W A = 


W A = 


Hv 


wu = 

m 


inertial angular velocity of body X, 

angular velocity of body X relative to body J(X), 

inertial angular velocity of momentum wheel m, 

angular velocity of momentum wheel m relative to body MO(m), and 


H m = angular momentum of wheel m about its own center-of-mass, relative to body 
MO(m). 

From definitions provided in virtually all texts on rigid-body dynamics, it follows that: 
a. Body X, a point mass, rigid body, or gyrostat 


II 

3 

>* 

-4* 

f 

(20) 


(21) 

= m K\K‘ 
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b. Body X, a point mass 


Li.x = T'i.x x m x‘ )r i,A 


L i,x ~ T'i.x x m x 1r i.x‘ 


c. Body X, a rigid body 


Li , a = *x ‘ 3 a + ^i.x x m x ^i.x 

Li.* " \ * "x + 3 x x *x * 3 x + ?,.x x m x\x- 
d. Body X, a gyrostat 

V * 

+ L*i H„ 


L,.* = *x * w x + T'l.x X m xT r u 


m: 

MO(m)= X 


v* 


L,.x = *x ’ 3 x + "a X ♦a * 3 x + ?i.x x m x ^i.x + “ H * 

MO(m)=A 


e. Momentum wheel m 


H m " *w * w w 

m nt 


H. = H m + I* -S, 

m 


H m = H m + I w * «x + "x X *w_ ' 


u, 


(22) 

(23) 

(24) 

(25) 

(26) 


(27) 


(28) 


(29) 


(30) 


It has been assumed that the composite N-body system has at most six degrees of freedom 
relative to an inertial reference, contiguous bodies have at most three degrees of relative 
freedom, and momentum wheels have only one degree of relative freedom. Thus, there 
exists at most 3(N + 1) + M degrees of freedom for the N-body system, and its motion 
can be completely defined by exactly N + 1 + M coupled-vector equations. 

These coupled-vector equations are given by the following: 
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a. For keS R 


This condition implies that body k is either a rigid body or a gyrostat and that the 
nest k-1 has at most three degrees of relative rotational freedom. The interaction 
force between bodies k and J(k), which constrains relative translational motion, is 
given from Equation (IS) by 



♦il-Ea.. 


Ff). 


(31) 


The rotation equation for all nests k-1 having keS R is obtained by substituting 
Equations (21), (27), and (31) into Equation (16): 




A • 5 a + WX «\7 1>x + 2 H„ 

in . 

— M°( m )=X 

- ♦ E K x • “x ff * ?r>l 


(32) 


b. For keS L 

This condition implies that the nest k-1 contains only the point mass labeled body k, 
and that the nest k-1 has at most three translational degrees of relative freedom. The 
translational equation for the nest k-1 is from Equations (IS) and (21): 


m, 


,7,., - Fj., ♦ F ; , ♦ Ff>- 


(33) 


c. The composite system is free to translate relative to the inertial reference. The 
translation equation for the composite system is from Equations (IS) and (21 ): 


m i = pH + pc + 2mJ F (e) . 
AeSp 0 0 tr 0 * 


(34) 


d. Momentum wheels are defined to have one degree of relative rotational freedom. 
The relative momentum equation for momentum wheel m contained in body X, 
derived from Equations (1 7) and (30), is given by 

“« + V ' 3 a = x *w * “a + CL m + ( 35 ) 

mm m 

These equations of motion must now be put into a computationally efficient form. 
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Accordingly, the inertial acceleration of the center-of-mass of body X can, by making use 
of Equations (2), (3), ( 1 0), and ( 1 1 ), be written as 






= fi , + 

* 1-0 


|3 I#) X D, . 3 X (3 X p,)| 

i* 1 


"x * “A + «x x K X 


XeS„ 


(36) 


“x + ^(x) x “x + 2w J(X) X ft X + ^JCX) X ( W J(X) X “x) XeS L» 


where the open dot implies differentiation with respect to the reference frame X fixed at 
hinge point X-l. 

Considerable symmetry in the final form of the equations of motion is obtainable if the 
inertial angular acceleration vectors are expressed in terms of the relative angular 
acceleration vectors. From the definitions of inertial and relative angular velocities 


CO. 


us 


& XeS D , 


(37) 


o,x-i 


and hence 


V' 


^ < 0 . 
fcS o,x-i 


(38) 


Substitution of Equation (38) into (36) along with a rearrangement of terms leads to 


7..X * 0, + 


V 


XeS„ 


us 


O.X-l 
i* X 


' v vy “► . — ► 

w ; X 7i-i,x + a x 


XeS, 
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w j(i) x ^ w j(o * 

i 

i# i 



«X X (w x X Q x ) 


XeS, 


(39) 


^ w j(\) X ® x + w j( X ) X (w J(X ) X a x ) XeS L - 
It is also convenient to make use of Equations (2) and (28) to write 


H m ~ V; w w m + w \ X H m , 


(40) 


where MO(m) = X and, again, the open dot implies differentiation with respect to the body 
X fixed reference frame. 

Substitution of the vector identities given by Equations (38), (39), and (40) into the 
system equations of motion, along with a rearrangement of terms yields: 

a. Rotation equation for the nest k-1 , keS R 


\eS 


k-1 


?k-i.x X m x 0, 


y z 


*W' w w 
m: mm 

MO(m)=X 


ieS 


0,k-l 

i#k 


y i y 

<I> x • co. + w 


\eS k{ 

L 




k-1 


Vu X m x (w | X Y,., iX ) 


ieS 


k-1 


4> x • Wj + ^ 7 k ., 


\ X m x (w f X 7,., x ) ieS R 


x*s, 




ii 


\u x m A 


ieS, 
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XfS 


k- 1 




• w x + / 


w x x '•'x w x 


m: 

MO(m)*X 


s » x »„ 


AfS 


k-l 


^k-I.X X m A 


kS 


°>}Q) X ^J(i) X 


O.X-l 

i*l 


^ x <"» x S x> 


2 W J(X) X Q X + W J(X) X ^ W J(\) x *V 


XeS„ 


XfS 


Lj 


^k-i + ^-i + 2 [”y k- 1 ,x x + • 


( 41 ) 


b. Translation equation for the nest X- 1, XeS L 


m. 


?i + x ?m.x + 

“ s o,x-i 

i#X 


-nu 


ieS 


Wjm X ( W J(|) X $ { ) + w,„v X (w J(x) X a x ) + 2w, w X 


o,x-i 


L 

+ h, + Fx H , + %>• 


ioo x kj 


( 42 ) 
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c. Translation equation for composite system 


m x@i x ?m.x) ieS R 


^S 0 " * KS 0 


m x a \ 


XL, X 


AeS ft ~ kS 


w j(i) x ( w ja) x + 


“x x <*x x *x> 


2w J(X) X a X + W J(X) X < W J(X) X a x) XeS L 




d. Momentum equation for wheel m:MO(m) - X 


w„. +£ 


V * «| " A x <V. ' 5 x + H m) + CL m + *w • <44) 


ff _ TT _■ . a ^ A O ▼▼ „ AM 111 \ ’ 

in m ieS A . , m m m 

0 ,A*i 

It is desirable to condense these vector equations eventually into a concise matrix format. 
To do this, however, several vector operations must be replaced by equivalent vector- 
dyadic, scalar-product operations. In particular, one can make use of the following: 

a. Multiplication of a scalar and a vector 

* m x 1 • t t , (45) 


m x“x * m x 1 • «x> 


where 


Is unit dyad. 


b. Vector product of two vectors 


w i x T'm.x " r n.x ‘ w , * 


j 



where 


T'k-l.X X "$1 " **Vl,X * 0, 

0 O 0 0 

(48) 

T'k-l.X X ®X * '^k-l.X * a X> 

(49) 


(50) 


and 3P is the tensor operator which transforms vectors into skew-symmetric tensors 
of rank two (dyads), if is an arbitrary vector with components | v,, v 3 , v 3 1 , then 


S>$) 


0 V 3 - V 2 

* V 3 0 V l 
V 2 -V, 0 


c. Vector triple product 

7k-l,X X "*X <«| X 7,1.X) * G k-1.M ‘ » 


where 


G k-i,n = -x K7 k .,.x ' T'n.x) 1 ' 7 n,x ^k-i.x) • 


(51) 


(52) 


(53) 


To establish symmetry conditions, it must be noted that the components of the 
pseudo-inertia tensors G£., and G£, k I >when written in the same coordinate system, 
obey the following matrix transpose relationship: 


l°k-l.i-|l “ l°M.kJ * 

The vector-dyadic form of the equations of motion is obtainable by substitution of 
Equations (45) through (52) into Equations (41) through (44). That is: 

a. Rotation equation for the nest k-1, keS R 


(54) 


“ m x r k .j x • ^ 

* 8 k-i ^ 1 * 8 o,k-i 


l*k 


XeS 

*S, 


k-1 


XcS 


k-1 


"k-U-l 


CO, 


i 

i 
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ieS 


k- 1 


£ ,.E 


**1.1 

A 


q* 

K ”, k,M 


S, KS n 


-E 


X£S ] 


nu r, 


1*1 


X*k-l,X M X 


E E 


lw • 

XeS. , m: m "m 


k-1 


ieS. 


MO(m)=X 


AfS 


k*l 
VS L 


+ Zj 


“x X \ ■ "x + ^ W X X H m 

m : 

MO(m)=X 


-E 


X<S 


kl 


^k-I.X X m X 


i«S 


W J(i) X ( W J(i) X 


0.X-1 

1*1 


*x x <*x x *x> 


XeS„ 


2 w j(\) x “x + W J(X) x ( W J(X) x a x) 

L 

♦ ♦ JJ, + E k.,., x f« * <•>] 

k-1 

b. Translation equation for the nest X-l, XeS L 


m, 


1 ■ ^ + us r u.x • "l + 1 • ? x 
ieS o,x*i 

i#X 


- nu 


(55) 


US Wj ® X X ^ + W, < X > X ^(M X + 2 W I(X) x 

""o.XI 

l#l 


l(X) " «x 


+ j$* + F h + 

X-I r X l T r x • 
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c. Translation equation for composite system 


^ m. 1 • <f. + ^ ^ 

X*S„ * ieS. 


X«S 


1-1 


m x r M,x’ W l ifS R 


*e S 


1-1 


"»X 1 ‘ fl X ifiS L 


m. 




Z u 

kS 0,k-l 

1*1 


J(0 


x (3 J(i) x ? ( ) + 


*x X ( 3 x x “x) 


XeS, 


20 J J(X) X tt X + W J(X) x ^ W J(X) x a x^ * cS l 


+ F* + F” + ^ F! 


F<*>. 




d. Momentum equation for wheel m:MO(m) = X 

4 y * 

*W m * S W m + . " «W m • W i 


(57) 


U X X 0w m * W X + H m) + CL m + * 


(58) 


EQUATIONS OF MOTION (MATRIX FORM) 

In the previous section, the equations of motion of the coupled N-body system have been 
derived in vector-dyadic form from the principles of linear and angular momentum. A 
cursory examination oi tquations (55) through (58) is sufficient to note that terms 
associated with the same gyroscopic effects have been grouped together. In putting the 
equations into a matrix format it is convenient to make use of partitioned matrices so as 
to retain this separation of effects. 

Accordingly, one may define: 

{« } = (N + 1) X 1 column matrix of relative acceleration vectors. The vector 
element in row k is 


a. Body k a rigid body 


w,. 


Wl 
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b. Body k a point mass keS L 


w. 


c. Total system k * N + 1 


‘k- 


(60) 


w k 


(61) 


| w w |= M X 1 column matrix of momentum-wheel relative angular acceleration 


vectors. The vector element in row m is 


co 


w 




(62) 


M- (N + 1) X (N + 1) symmetric matrix of pseudo-inertia tensors. The tensor 
element in row k, column i (i > k) is from Equations (55), (56), and (57), given by 

a. keS R ieS k ., ieS R (Nest i-1 is a nest of at least one rigid body within the nest k-1). 


k,i 


*X + ^ ® k- 1 ,1- 1 • 


Xf s 


M 

\es. 


XeS 


(63) 
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For k = i, X kJ is the inertia tensor of the nest k-1 about the hinge point k-1 . For 
k ¥> i, ii defines an inertia cross-coupling te- -or between nests k-1 and i-1. 


b. keS R ieS k , ieS L (Nest i-1 is a nest of one point mass within the nest k-1.) 


X|c,| “ ‘ m i ^k- 1 ,i • 


(64) 


is the tensor form of the mass moment of the nest i-1 relative to the hinge point 
k-1. 


c. keS Ri i#S k , (Nest i-1 is not contained within the nest k-1.) 


*k.i = °- 


d. keS L i * k (Nest i-1 is identical to nest k-1 and contains one point mass.) 


x k | -111,1. 


X k(i is the tensor form of the mass of the nest i-1. 


(65) 


( 66 ) 




22 



1 


e. keS L i^S k , (Nest i-1 is not contained within the nest k-1.) 

f. keS_ i * N + 1 (Nest k-1 has at least one rigid body.) 






X kJ is the tensor form of the mass moment of the nest k-1 relative to the hinge point 


g. keS L , i * N + 1 (Nest k-1 has one point mass.) 


X k,l “ m k 1 • 


X M is the tensor form of the mass of the nest k-1. 


h. k“N + l,i-N + l 


m. 1 , 


keS„ K 


X kJ is the tensor form of the total system mass. 

It should be noted that when the tensors defined above are all expressed relative 
to a common frame of reference, the matrix (X] is symmetric. This is easily 
proven by application of Equations (SO), (51), and (54). 

L''J = (N + 1) X M, rectangular matrix of momentum wheel inertia tensors. The 
tensor element in row k, column m is from Equations (55) and (56) given by 

a. MCXm) eS k l (The rigid body in which momentum wheel m is imbedded is within 
the nest k-1.) 

,c k.m * I w m * (71) 

b. MO(m) *S k ., or k * N+l (Momentum wheel m is not within the nest k-1 .) 


■u - o 


U*J T * M X (N + 1), rectangular matrix, the transpose of (j c J , 

M • M X M, square matrix of momentum wheel inertia tensors. From 

Equation (58), the tensor element in row m, column n is given by 


r •- 

i 



a. m * n 


* ‘ w m * (73) 

b. m # n 

^m,n * ® * (74) 

{»?'} s (N + 1 ) ,' 1, column matrix of forces and torques associated with 
centripetal and Coriolis acceleration effects. Note that the force 
associated with the mass of body X and its combined centripetal and Coriolis 
acceleration can be expressed as 


c x " “a 


~ x (3„ x li> ♦ | 

""O.A-l 
I# I 


“x X (3 X X a x ) 


2w j(X) X a \ + W J(X) X ^ W J(X) X a \ ) 


XeS R 

XeS L 


(75) 

From Equations (55), (56), and (57) the vector element q? ir. iow k is 
a. keS R (The nest k-1 h?.s at least one rigid body.) 





(76) 


tj' is the resultant torque at hinge point k-1 due to the centripetal and Coriolis 
acceleration of each body in the nest k-1. 

b. keS L (The nest k-1 contains one point mass.) 

- -C, • (77) 


c. k * N + 1 



XfS 


0 


(78) 


r?J, +| is the resultant force acting on the composite system due to the centripetal 
and Coriolis acceleration of each body. 


,-,<S tiw 
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(N + 1) X 1, column matrix ot torques associated with the inertial 
angular velocity of the body fixed reference frames and the inertial 
angular momentum of the rigid bodies and gyrostats about their respective 
centers-of-mass. Note that the inertial angular momentum of body X 
about its own center-of-mass is given by 




♦» 


* w. ♦ ^ H . 


m: 

MO(m)*X 


From Equation (55) the vector element tj}, in row k is 


(79) 


a. keS R (Nest k-1 contains at least one rigid body.) 


n 


i 

k 


^ * i k,\ • 

* s k-> 


(80) 


Tj l k is the rate of change of angular momentum of all rigid bodies and gyrostats in 
the nest k-1 about their own cente.s-of-mass due to the inertia) angular velocity 
of their respective body h/ed-reference frames. 


b. keS L or k * N + 1 


n 


i 

k 


0. 



(81) 



M X 1, column matrix of torques associated with the inertial angular 
velocity of the gyrostats fixed reference frames and the inertial angular 
momentum of each momentum wheel about its respective center-of-mass. 
Note that the inertial angular momentum of momentum wheel m imbedded 
in body X about its own center-of-mass is given by 



.w 


m 


'"w 





(82) 


From Equation (58) the vector element in row m is 


< ■ -X * • <83) 

tj * is the rate of change of angular momentum of momentum wheel m about its own center- 
of-mass due to the inertial angular velocity of the reference frame fixed in the gyrostat 
in which it is imbedded. 


S 


y 
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0 C * = (N + 1) X 1, column matrix of the forces and torques of constraint, 

4 1 = (N + 1) X 1, column matrix ol the forces and torques associated with 

' mechanisms existing between contiguous bodies, and 

f E ^ 

|0 1 J = (N + 1) X 1, column matrix of the forces and torques associated with 
causes external to the N-body system . 


The vector element in row k of 


{♦'■} * {*"■! * !♦'■} 


a. keS B (Nest k-1 contains at least one rigid body.) 


V * * *1' - K. 




This is the resultant torque external to the nest k-1 about the hinge point k-1. 

b. keS L (Nest k-1 contains one point mass.) 

*5' * <' * <■ - n„ * * ?<•>. 


This is the resultant force external to the nest k-1. 


c. k = N + 1 


K'**? 1 + = K + ‘•o 


This is the resultant force external to the composite system. 


|0 C2 1 = M X 1, column matrix of the constraint torques acting on the momentum 
wl.eels, and 

U Hj | = M X 1, column matrix of the torques associated with mechanisms existing 

' ' I_ _ it ...L t. A iU a ...U ' AM 


be.-veen the momentum wheels and the bodies in which they are 
imbedded. 


The vector dement in row m of 


(* c ’l * M 


V 



is 


*3 * ♦£• • \ ♦ CL. . (| 

This is the resultant torque external to and acting on the momentum wheel m. 

Making use of the above notation, the simultaneous vector-dyadic equations of motion 
given by Equations (55) through (58) may be expressed in partitioned matrix form as 



To solve these equations numerically, they must be put in the form of a set of 
simultaneous, scalar-differential equations. The difficulty with this step is in solving for 
or eliminating the unknown forces and torques of constraint. 

In theory, if the N-body system has a total of N T degrees of freedom and N L holonomic 
conditions of constraint, then it is possible to generate a set of 

N f = n t - N l (89) 

generalized coordinate equations which are independent and completely define the system 
dynamics. 

In Reference 2, Hooker and Margulies present a method for constraint elimination. The 
method essentially derives expressions for the constraint torques in terms of system 
parameters and then substitutes them into the equations of motion. From a computational 
viewpoint, this method is cumbersome since the order of the system of equations 
remains unchanged and several additional matrix inversions are required. 

In Reference 3, Hooker presents another method for constraint elimination which ’.s 
well suited for digital computation. The method defines a set of N F free coordinate 
vectors which span the N F dimensional vector space in which motion is possible. By 
application of orthogonality conditions, forces and torques of constraint can be 
eliminated and a set of N F simultaneous scalar equations derived which completely 
define the system’s motion in the N F dimensional vector space. 

To accomplish this, Hooker notes that at every hinge point it is possible to define a triad 
of linearly-independent free and locked coordinate vectors. Physically, if one imagines 
that contiguous rigid bodies are connected by either a zero-, one-, two-, or three-axis 
gimbal, the free coordinate vectors correspond to unit vectors fixed along the gimbal 
axes, and the locked coordinate vectors correspond to unit vectors fixed along axes about 
which motion is totally constrained. 


27 



Since the equations of relative motion are vector equations, at each hinge point the 
relative velocity vector is expressible in terms of a linear combination of the free 
coordinate vectors while the constraint torque is expressible in terms of a linear 
combination of the locked coordinate vectors. These vectors are orthogonal to each other. 
By forming the appropriate vector scalar products, the scalar equations which define 
the components of motion about or along each free coordinate vector can be obtained 
and the need for evaluating the 'Tmstraint torques circumvented. 

To apply Hooker’s technique, several kinematic related definitions must be introduced. 


KINEMATICS OF THE N-BODY SYSTEM 

The subject of kinematics to which this section addresses itself is concerned exclusively 
with determination of the relative orientation and rate of each body in the N-body system. 
The causes to which relative motion may be ascribed are not of interest. User convenience 
and computational efficiency have been the primary factors used in determining which 
kinematic methods should be applied. 

It is the author’s opinion that there is no one kinematic technique which is best for all 
possible problems. Accordingly, the program N-BOD has been programmed to give the 
user a limited selection of kinematic options. Unless directed otherwise, direction cosine 
methods are applied. No attempt is made to artificially orthonormalize the computed 
transformation matrices. At the user’s option, an algebraic quaternion method can be 
used. This is simply a generalization of standard Euler angle techniques. Selective use of 
both direction cosine and algebraic quaternion methods should permit the user to 
circumvent the common kinematic problems of orthogonalization and gimbal lock. 

At the inertial origin and at every hinge point of the N-body system a reference frame 
has been defined; 

X® X® X® | = reference frame 0, fixed inertially at the inertial origin, and 

X* Xj X* | = reference frame X, fixed at hinge point X-l of body X; 

if body X is a point mass, the coordinate axes are respectively 
parallel to those of reference frame J(X). 

It is desirable to provide the user with the ability to specify an arbitrary “at rest,” or 
“nominal zero internal stress” orientation for each reference frame. Accordingly let 



transformation matrix which takes vectors from body J(X) fixed 
coordinates into body X fixed coordinates when the system is at rest 
in the nominal zero internal stress state. 


To completely define the at-rest state of the N-body system, translation aH momentum 
wheel orientation conditions also must be defined. 
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a. The entire N-body system is ftee to translate relative to the inertial origin; 

if, = vector from the inertial origin to the origin of reference frame i 
(the center-of*mass of body 1). 

b. Point masses may be given an attest position relative to their respective 
contiguous bodies; 

•*> 

N a x = vector from hinge point X-l to the at-rest position of the poir 
massX. 

c. The rotation angle through which a momentum wheel has rotated can be 
important for special cases. Rather than introduce a wheel-fixed reference frame, 
it is defined that the at-rest angle of the wheel, relative to the body in which it is 
imbedded, is zero. There is no loss of generality here since all momentum wheels 
are by definition symmetric. 

In order to eliminate the forces and torques of constraint in the equations of motion, a 
triad of unit coordinate vectors is defined at every point in the system about or from 
which relative motion is measured. 

Let 


N 0 = total number of coordinate vector triads. For a system of N coupled 
bodies and M momentum wheels 

N 0 = N + 1 + M. (90) 

N p = total number of free coordinate vectors; these span the N p dimensional 
vector space in which motion is possible. 

total number of locked coordinate vectors; these span the N L dimensional 
vector space in which motion is totally constrained. 

■ • > Qn f } - set unit free coordinate vectors which span the N F 
dimensional vector space in which motion is possible. 

• • > Pn l | = set ot unit locked coordinate vectors which span the N L 
dimensional vector space in which motion is totally constrained. 

. .,0 Nf J * set of scalar parameters which at any instant of time 
define the relative displacement about or along each respective free 
coordinate vector (analogous to Euler angles). 

. . , 0 Np | = set of scalar parameters which define the time rate of 
change of the respective displacement parameters (analogous to 
Euler angle rates). 


(?,A- 

(».».,■ 

K «». ■ 
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To uniquely define each of the free and locked coordinate vectors, several definitions 
must be made. In principle the vectors may be randomly labeled; however, for 
computational purposes it is advantageous to choose a particular numbering st ance. 

Each triad of coordinate vectors may be assigned a label. The origin of the triad is the 
point from or about which relative motion is measured. That is: 

a. To measure the relative rotation of rigid body k, the triad k is defined to ha.. 
its origin at hinge point k-1. 

b. To measure the relative translation of point mass k, the triad k is defined to have 
its origin defined at .is at-rest position. 

c. To measure the relative translation of the center-of-mass ol oody 1 relative to the 
inertial reference, the triad N + 1 is defined '<> have its origin at the inertial origin. 

d. To measure the relative rotation of momentum wheel m, the triad N + 1 + m is 
defined to have its origin at the center-of-mass of mor lentum wheel m. 

The mix of free and locked coordinate vectors in each triad is defined by the integer 
function P (X) where 

P(X) = total number of locked coordinate vectors in the triad X. 

The free and locked coordinate vectors which make up triad X are: 

a. P(X) = 0 Three degrees of freedom 

{ q j ,q j+» ,q j+2 }> 

b. P(X) = 1 Two degrees of freedom 

{ q i-V.’P.}’ 

c. P(X) = 2 One degree of freedom 

q j’ Pi’ Pi+l J * an d 

d. P(X) = 3 Zero degrees of freedom 

where for any X, X = 1, 2, . . . , N„ 

j = 1 + 


s 


(3 - P(fi» , 


(91) 
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i - 1 + z* m> 


(92) 


and 

= 0 . (93) 

If body X is a point mass, the free and locked coordinate vectors of the triad X are 
mutually orthogonal and fixed with respect to the reference frame X. 

If body X is a rigid body, then by definition it is connected at hinge point X-l to body 
J(X) by either a zero-, one-, two-, or three-axis gimbal. The numerical ordering of the free 
coordinate vectors is defined by the following convention: 

a. For a one-, two-, or three-axis gimbal, the first rotation is about qj through the 
angle 6 .. The free vector qj is fixed with respect to the reference frame J(X). 

b. For a two- or three-axis gimbal, the second rotation is about <f. +1 through the 
angle 0. +1 . For a two-axis gimbal, the free vector qj +1 is fixed with respect to 
reference frame X. For a three-axis gimbal, the free vector q^ , is defined by the 
vector cross-product 




(94) 


c. For a three-axis gimDal, the third rotation is about q^ 2 through the angle 0 . +2 . 
The free vector qj +2 is fixed with respect to the reference frame X. 


If body X is a rigid body, the numerical ordering of the locked coordinate vectors at the 
hinge point is arbitrary. With the exception of the case of a two-axis gimbal, locked 
coordinate vectors are fixed in reference frame X. For a two-axis gimbal, the locked 
vector pj must be orthogonal to both free vectors qj and qj +i . Thus, it is defined by the 
vector cross-product 


x 

9] * V, 

If X = N + 1, the free and locked coordinate vectors are fixed with respect to the inertial 
reference frame. They are defined by convention to be aligned respectively with the 
X®, X®, and X° coordinate axes. 
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If X = N + 1 + m, the free vector is aligned with the spin axis of momentum wheel m which 
is fixed in reference frame MO(m). The two locked coordinate vectors are fixed normal 
to each other and to the spin axis in a wheel fixed-reference frame. 

INITIAL ORIENTATION OF THE N-BODY SYSTEM 

In defining the relative orientation of the bodies in an N-body system at time zero, it is 
undesirable to start from the hypothesis that the elements of the initial transformation 
matrices are given. If the response of a complex system under different initial conditions 
is to be studied, the computation of initial transformation matrices can be tedious. It is 
desirable to define a set of physically realizable, independent parameters which can be used 
to construct internally the initial conditions necessary for computation. 

The scalar parameters which define relative displacement and rate about and along the free 
coordinate vectors provide a desirable set of parameters to work with. 

Since relative motion takes place only about or along the free coordinate vectors, the 
initial orientation can be given by stipulating exactly N F independent relative displacement 
parameters, one associated with each free coordinate vector. These parameters are taken 
relative to the nominal zero stress position. 

Let 



transformation matrix which takes vectors from body J(X) fixed 
coordinates to body X fixed coordinates at time zero. 


To compute all initial transformation matrices [° ^ = 1, 2, . . . , N, it is not 

possible to assume that all gimbal axes will be parallel to body fixed-coordinate axes. It 
is therefore most convenient to make use of the quaternion techniques reviewed briefly 
in the appendix to compute the initial transformation matrices. 


Making use of the matrix operator & and the quaternion operator 3t defined in the appendix, 
the initial value of the transformation matrix [ , X = 1 , 2, . . ., N is given by 


a. One-axis gimbal, reference frame Ij aligned with body X fixed axes 



(96) 


b. Two-axis gimbal, reference frame I 2 aligned with body X fixed axes 


K'VJ = '^ (> N f, i h\) 


r N ->~ l 

Lx J(X)J, 


(97) 
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c. Three-axis gimbal 

s '*(i N Si t ijfij iJ\J 

where 

J,, ((?„),,) 



= cos0 m+1 /2 + <? ^ j 

^m+1 


(100) 


= cos0 m+J /2 + j 

^m+2 


(101) 


tw, 


sinfl /2 


(98) 


(99) 


SYSTEM ORIENTATION VIA DIRECTION COSINES 

To integrate the equations of motion of the system it is convenient to refer all vectors and 
dyads to a common frame of reference. Depending upon application, this may be either the 
inertially fixed frame of reference or the coordinate axes fixed in body 1. Generality is 
retained since the system can always be relabeled so that body 1 does contain the com- 
puting frame of reference. 

Let 

[>1 = transformation matrix which takes vectors from the computing frame 
L of reference into body X fixed coordinates at time t; t > 0, X = 

1, 2, .... N. 

— ► 

To transform the vector R given relative to the computing frame coordinates into body X 
fixed coordinates apply 

{*}x "[*<]{»}«• (102) 

To transform the tensor R given relative to the computing frame coordinates into body X 
fixed coordinates apply 


[«ix - (»i. (A), 


where 




(103) 


(104) 
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If the computing frame is chosen to be the inertial-reference frame, then c * 0 and 




. TT 


keS 


lk ,y j(k)l f 


0,\-l 


( 105 ) 


where 



[k'^I(k)] J(X)’^J(J(\))] • • • [r^o] • 


(106) 


If the computing frame is chosen to be the body 1 fixed-coordinate frame, then c = 1 and 


- TT 

A 1 k*S„ * 


0,\-i 

k*l 


lk -:/ J(k)] • 


(107) 


Since the components of the N-direction cosine transformation matrices [ r / c ] are time 
varying, they must be continually updated. An expression is required to define the time 
rate of change of the components of [ x ./ c ] • These quantities can then be integrated 
and the time histories of the matrices defined. 

Let 

R = arbitrary vector fixed in body X, and 

c co A = angular velocity of the body X fixed-reference frame relative to the 

computing frame of reference. 

From vector differential calculus 


c dR x dR -> 

— = — + c w x XR » (108) 


where the presuperscript on the d/dt operator defines the reference frame in which 
differentiation is referenced. But R is fixed in body X, therefore 


c dR _ ■* 

— = A x R - 

dt c A 


(109) 
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T 


w 


column matrix of the components of R relative to the computing 
frame, and 

column matrix of the components of c w x relative to computing frame 


Recall that if 


(a) < ■ 


CO 

1 


( 110 ) 


then the matrix operator SP can be defined such that 

1M.)- 


SP\ 


0 W 3 -W 2 


-w 3 0 U, 

u> 2 -to, 0 


( 111 ) 


and vector Equation ( 1 09) can be written in matrix notation as 

PI. ~'(Ia|.) 1*1.. 

where the dot over the matrix implies differentiation of each matrix element with 
respect to time. 

From the definition of [ x .^' c ] note that 

{*}»* w {»},. 

Since R is fixed in body X 

{S)x - 0, 
and 

(«}« ■ W T {*}x ■ lAl T u-^i {*},. 


( 112 ) 


(113) 


(114) 


(115) 
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It follows from Equations (1 12) and (1 14), since the vector R is arbitrary and fixed in 
body X, that 

~ * jc^- ( 116 ) 

Take the transpose of both sides and premultiply these by [ x '^ c ] to obtain 

U./-J = //((a),). on) 

Note that since 

= W, (118) 

Equation (117) can also be expressed as 



Equation (117) defines nine differential equations, one for each direction cosine. To 
completely define the relative orientation of the N-body system, this implies that 9N 
direction cosine equations must be solved simultaneously with the system equations of 
motion. 

Due to the orthonormality of the direction cosine matrix, only six of the nine equations 
need be integrated. The solution to these equations can then be used to compute 
algebraically the three remaining quantities. 

SYSTEM ORIENTATION VIA QUATERNION METHODS 

The computation of transformation matrices may be accomplished by either integrating 
a set of kinematic differential equations, such as direction cosines, or algebraically by 
application of quaternion methods. Rather than set up a series of quaternion differential 
equations, as is commonly done, it is possible to make use of the free coordinate vectors 
and the respective angles of rotation to simply construct the required transformation 
matrices algebraically. It becomes a debatable question as to which technique is 
optimum for all problems. 

To apply quaternion methods, the same procedure used to obtain the initial orientation 
of the N-body system is used. 




That is, from 


and Equations (96) through (101) it follows that 
where 


'A 

one-axis gimbai I 1 = X 

A A 

two-axis gimbai l 2 = X 

A A A 

three-atis gimbai t 


(120) 


( 121 ) 


( 122 ) 


It should be noted that the above equations reduce to one, two, or three successive 
Euler-angle rotations when the free coordinate vectors are aligned with coordinate axes 
in the at-rest state. 


RELATIVE RATE AND DISPLACEMENT 


At any hinge point k-1, the relative angular or linear rate between the contiguous bodies 
J(k) and k can be expressed as 


u 


k 



0 23) 


where the summation is carried over all free coordinate vector indices defined at the 
origin of the triad k. 

If body k is either a point mass having one, two, or three degrees of relative translational 
freedom or a rigid body having one or two degrees of rotational freedom, the relative 
displacement parameter 0 m+ . along or about the free vector q > m+j is obtained from the 
equation 


^m + l 



(124) 
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where 


^m + l " ^m+l ' w k *“0,1,2. (125) 

If body k is a rigid body hinged to body J(k) by a three-axis gimbal, then since the three, 
free-coordinate vectors are not always mutually orthogonal, the relative rotation param- 
eters 0 m+j (i * 0, 1,2) must be derived as follows. 


Let 

-► 

9m 

* "k * 

“ 9 m 

m m 

+ c ^m+2 • 

(126) 



-► £ 

‘ W k 

* n m + . “ 

®m+l * 

(127) 


9 m +J * W k 

“ c6 m 

+ 9 m*i » 

(128) 

and 


9» 

• %,+J “ c - 


(129) 


From a simultaneous solution of the above 


m + l 



i * 0, 1, 2 , 


where 


and 


■ cn m + 2 


1 “ V 


i + l * ^m + l * 


m + 2 


1-C J 


ELIMINATION OF CONSTRAINT TORQUES 


(130) 


(131) 


(132) 

(133) 


The partitioned matrix form of the simultaneous vector-dyadic equat’jns of motion for 
the coupled N-body system are from Equation (88) given by 
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f * 


r * 
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To obtain a solution of this equation, the forces and torques of constraint, defined by 
the elements of the matrices j <j> 'J and {0 2 f, must be either analytically defined or 
the equations restructured so as to eliminate them. The latter approach, used by Hooker, 
is adopted here. 


To restructure the equations and eliminate the necessity for evaluating the forces and 
torques of constraint, the procedure outlined on page 27 in the section entitled Equations 
of Motion (Matrix Form) is followed. 

From Equation (123), 


w. 


m (a k- i 


0 q 
m ~ m 


k = 1, 2 N, N + 1 


( 134 ) 


where k = N + 1 implies the inertial origin. In a parallel manner 


'w 




in = 1, . . M , 


where the index j is defined by Equation (91 ). 
Differentiation yields 


and 



(0 ? + 0 cT 1 

m m "m 1 


O 



(135) 


(136) 


( 137 ) 


where the closed and open dots imply, as before, differentiation with respect to the 
inertially fixed and locally fixed reference frames, respectively. 


In order to put Equations ( 1 34) through (137) into the matrix format required for 
substitution into Equation (88), the following definitions are made. 

Lq J = (N + 1) X (N f - M) rectangular matrix of the free coordinate vectors. 
The vector element in row k, column m is 

a. k = 1.2, ...N 


^k.m 


if free vector m is defined between bodies J(k) and k 
if not . 


(138) 
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b. k = N + 1 


9k, m 


q m if free vector m is defined at the inertial origin 
0 if iot . 


U39) 


[h] = M X M square matrix of thr free coordinate vectors existing along 

momentum-wheel spin axes. The vectoi element in row m, column n is 


h 


m ,n 


{ ‘•Np -M+m 
0 


if m - n 
if m # i. . 


(140) 



(N F - M) X 1 column matrix of scalar rates about or along the free 
coordinate vectors existing between bodies and at the inertial reference. 
The scalar element in row m is 0 . 

m 


M X 1 column matrix of momentum-wheel relative rates. The scalar 
element in row m is ff w . 

in 


Making use of the matrix notation, one mav write 


and 



(141) 


(142) 


To eliminate the forces and torques of constraint, note that at every point from or about 
which relative motion is measured, the constraint vector is normal to all free coordinate 
vectors defined there. Thus it follows that the matrix vector-scalar product equation 


A 4-°. 
1 

T 



0 


c. 



— ( 

O 

1 


U 2 


0 


( 143 ) 
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can be used to eliminate the forces and torques of constraint from Equation (88) and 
restructure it into a set of simultaneous scalar equations. 

Direct application of Equations (141) through (143) in Equation (88) yields 


IvIvf'BHvHvivJItl 


Lo i hj L cT t P J Lo ! 0 J l «J 


..-t 


0 i 


V 

+ — + 


\L 

^0 


(144) 


Equation (144) defines the N F scalar equations which completely define the motion of the 
coupled N-body system. These equations have been programmed for numerical solution 
and form the basis for the digital computer program N-BOD. 


By intent, the only elements in Equation (144) which have not been extensively discussed 
are the elements of the matrices 




The elements of these matrices cannot in general be generalized. They define the forces 
and torques due to mechanisms existing between contiguous bodies, such as springs, 
dashpots, motors, and so forth and the forces and torques associated with effects external 
to the N-body system. The inclusion of several effects often encountered will be discussed 
in the following sections. 

FORCES AND TORQUES DUE TO NONGYROSCOPIC EFFECTS 

Most problems of practical interest involve the inclusion of one or more nongyroscopic 
effects into the equations of motion. For example: 

a. Contiguous bodies may be viscoelastically coupled at their respective hinge 
points. 
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b. Control systems may exist which, based upon a given set of control laws, 
activate motors that drive contiguous bodies relative to each other or alter the 
angular rates of momentum wheels. 

c. Thrusters controlled by an active control system may exist at various points in 
the N-body system. 

d. The system may be disturbed by one or more forms of environmental loading 
such as that due to gravity, gravity gradient, thermal, solar pressure, 
aerodynamics, etc. 

It is a relatively simple matter to include any or all of these effects in the matrices 

MW 1 ) 


defined in Equation (144). 

ELASTIC COUPLING OF CONTIGUOUS BODIES 

Springs can be placed between any pair of contiguous bodies. These restrict relative 
rotational motion for rigid bodies and relative translational motion for point masses. The 
springs may be either linear or nonlinear; however, for simplicity, this discussion 
considers only linear springs which restrict motion about or along free coordinate axes. 

Consider the nest k - 1. Note that at every hinge point between bodies contained within 
the nest k - 1 spring forces and torques appear in equal and opposite pairs and thus 
exactly cancel each other. At hinge point k - 1 , however, the reactive spring torque 
(force) is external to the nest k - 1 and must be treated as an external disturbance. 

Assume that <f m is a free coordinate vector defined at hinge point k - 1 and that the 
motion about or along is restrained by a linear spring. 

Let 


K m = spring constant of the linear spring which restrains motion about 
or along the free coordinate vector ; 


then 

- K m 0 m cf m = spring torque (force) associated with the relative angular 

(translational) displacement 8 m of the contiguous bodies k and 
J(k) about (or along) the free coordinate vector <f m . 


It follows that the net reactive spring torque (force) at hinge point k - 1 is 


-E 

m@k- 1 


K m e q„ 


k = 1,2, . . . ,N + 1, 


045) 
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where the summation extends over all free coordinate indices m defined at hinge point 
k - 1 of body k, for k = N + i the free coordinate indices defined at the inertial origin 
are used. 

This vector quantity must be added to the row of j^* 1 ] • 

Similarly, if the momentum wheel n is elastically coupled to the body in which it is 
imbedded, the spring torque 


K m 0 w m q m’ 


(146) 


where 


m = N F - M + n, (J4 , 

must be added to the n ,h row of j/* 2 } . 

DISSIPATIVE COUPLING OF CONTIGUOUS BODIES 

Dampers can be placed between any pair of contiguous bodies to retard relative motion. 
For simplicity, only linear viscous-type damping mechanisms are considered. Nonlinear 
dissipative devices such as hysterisis or Coulomb friction dampers can be incorporated 
in the formalism; however, their modeling can become quite complex. 

Consider the nest k - 1. As in the case of springs, the only reactive damper force or 
torque acting external to the nest is the one defined at hinge point k - 1. 

Assume that tf m is a free coordinate vector defined at hinge point k - 1 and that the 
relative motion about or along cf m is retarded by a linear viscous damper. 

Let 


D m = damping coefficient of the linear viscous damper which retards 
motion about or along the free coordinate vector q" n> ; 

then 

- D m 6 m q m = damping torque (force) associated with relative angular (translational) 
velocity 9 m of the contiguous bodies k and J(k) about (along) the 
free coordinate vector ql . 

*m 

It follows that the net damping torque (force) at hinge point k - 1 is 

k-1.2 N+l. (148) 


•£ 

m@k- 1 




i 
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This vector quantity must be added to the k f * row of . 

Similarly, if the momentum wheel n is dissipatively coupled to the body in which it is 
imbedded, the damping torque 


- D_ 6 


m w W_ q m> 


049) 


where 


m = N F -M + n, 


must be added to the n ,h 


row of 



MOTOR COUPLING OF CONTIGUOUS BODIES 


(150) 


Motors may be used to actively control the relative orientations ot contiguous bodies; 
they may also be used to control the relative rates of momentum wheels. In most 
practical problems of interest, the motor torques are defined as the outputs of an active 
control system; the control system makes use of various system state variables and control 
laws to define the appropriate motor torques required to achieve a predefinable 
objective. 

Let cf m be the free coordinate vector between bodies J(k) and k about which a motor 
exists. Furthermore, let 

CL m = scalar magnitude of the motor torque to be applied by the motor 
about the free coordinate vector q ; 

‘HI 7 


then 

CL m Tm = motor torque applied about free vector q“ m to body k. 
It follows that the net motor torque at hinge point k-1 is 


L 

m@k- 1 


CL m9m 


k = 1,2, . . . ,N + 1 . 


(151) 


This vector quantity must be added to the Y* h row of |0 Hl J • 

Similarly, if the momentum wheel n is coupled by a motor to the body in which it is 
imbedded, the motor torque 


CL m % , 


(152) 
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where 


m * N F - M + n, 


must be added to the n th row of 



(153) 


LOCALLY APPLIED FORCES 

The type of locally applied forces which most often occur in satellite simulation are those 
attributable to gas jet firings. The location and force of these gas jets are definable. Let 

R x j = radius vector from the center-of-mass of body X to gas jet J which is 
located on body X. 

F^> = force associated with gas j~t J of body X when fired. 

It is perfectly admissible for the firing of these jets to be governed by a control law which 
is a function of the relative attitude motion of the system. It is also admissible for the 
force to build up as some definable function of time. 

From Equations (32), (33), and (34) it can be seen that the thruster acting on body X 
produces a force which is external to every nest of bodies containing the body X. The 
vector quantity which must be added to row k of \<1> ' } is 


0 

ifX*S k ., 

p( e ) 

r \,J 

if XeS k , and keS ] 


if XeS k and keS, 

p(<?) 

r \,J 

if k = N + 1 


(154) 


where k = 1, 2, . . . , N + 1 , 

DISTRIBUTED FORCE FIELD 

Environmental loading due to gravity, aerodynamics, solar pressure, and so forth produces 
a force field which is distributed over the N-body system. 

From Equations ( 1 8) and ( 1 9) let 

F* e) = resultant force acting on body X due to the distributed force 
field, and 



resultant torque acting on body X due to the distributed force field. 



Again, it follows from Equations (32), (33), and (34) that the vector quantity which must 
be added to row k of {0 E, f is 


c(e) 

r \ 


X = k, keS L 


^ j [>,.x x F l e> + keS R 


(155) 



p( e ) 

r x 


k = N + 1 . 


For the particular example of an earth based system subject to a uniform gravitational 
force field, the gravity force acting on body X is 


where 



(156) 


g = acceleration of gravity, 

0 o = unit vector directed from the earth’s center to the center-of-mass of 
the composite N-body system. 

Furthermore, since the gravitation force is distributed uniformly through the entire 
volume of each body, 


-*<e) 

<P\ = 


0 . 


(157) 


In this situation, /?, is measured from the inertial origin which does not necessarily have 

to be at the earth’s center. @ 0 is simply a unit vector, which can be expressed in any 
coordinate system that is constant in the inertial reference frame. 


ORBIT DEFINITION 

To incorporate into the formalism the ability to study gravity-gradient effects or to 
include earth, sun, or star sensors in an active on-board control system simulation, a 
rudimentary definition of the orbit must be made. 
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Assume that the satellite is in an elliptic orbit around a spherical earth (see figure 3). 
To define the orbit, the following quantities must be given: 

a = semi-major axis of elliptic orbit, 

e = orbit eccentricity, and 

T # = time of perihelion passage. 

From these quantities and elementary orbital mechanics, one may compute 
P 0 = distance from earth center to composite system center-of-mass, 

v = true anomaly. 



Let 

G e = earth’s gravitational constant; 
then the orbit period P e is 

P e 

the mean motion rj e is 
and the mean anomaly M e is 

K 



(158) 

(159) 


( 160 ) 
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From Kepler’s law, the eccentric anomaly E # is given by 


M, = E e - e sin E # . 

It has been shown in various texts that 


0 O = a (1 - ecosE # ) * 


a (1 - e») 
1 + e cos v 


( 161 ) 


(162) 


and 


tan 


v 

1 




(163) 


or 


cosv = 


cos E e - e 
1 - e cos E e ’ 


(164) 


and 


sinv = 


•J\ - e 2 sin E e 
1 - e cos E. 


(165) 


The problem is to obtain E. as a function of M„. By the method of Lagrange, Moulten 

€ v 

(Reference 8) shows that 


e* 

E = M + e sin M + — sin 2 M. 

c c c 2 " 

+ _£ (3 J sin 3 M e - 3 sin M e ) 

3 ! 2 2 

+ (4 3 sin 4 M - 4 • 2 3 sin 2 M„) 

4! 2 3 

e 5 

+ (5 4 sin 5 M, - 5 • 3 4 sin 3 M # + 10 sin M # ) 

5 ! 2 4 

+ (6 s sin 6 M # - 6 • 4* sin 4 M # + 1 5 • 2* sin 2 M # ) + . . . , 

6 ! 2 s 


( 166 ) 
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The expression converges rapidly for small values of eccentricity and can be suitably 
truncated in application. 


In any reference frame fixed in time, relative to an earth fixed-reference frame, the vector 
from the earth’s center to the composite center-of-mass of the N-body system may be 
written as 


t = ft 


o » 


( 167 ) 


A 

where the magnitude of 0 O is defined above and the components of the unit vector p 0 
are trigonometric functions of v, the true anomaly. 


GRAVITATIONAL FORCE FIELD FOR ORBITING SPACECRAFT 

For large orbiting spacecraft, it is improper to assume that the gravitational force is 
uniformly distributed through the entire structure. For a significant class of problems, 
gravity-gradient effects can influence the attitude dynamics of the spacecraft. 

— ►f. 

From Newton’s law of gravitation, the gravitational force F£, acting upon the mass 
element m. . of the body X in the N-body model of the spacecraft is 

A4 



= - G e m A,i 


'I,X,i 




X,i 


* 1,X, 


,) 3/2 


( 168 ) 


where 


G e * earth’s gravitational constant, 

= vector from the earth’s center-of-mass to the mass element m XJ . 

Let 

jf 0 = vector from the earth’s center-of-mass to the composite N-body system’s 
center-of-mass, 

* vector from composite N-body system’s center-of-mass to the center of 
mass of body X, and 

R XJ * vector from the center-of-mass of body X to the mass element m x . . 

Then 




( 169 ) 


49 


y 


V 



I 


t 


By substitution of a truncated binomial-series expansion for the denominator of 
Equation (168) it follows that 


_> r • .(\ + R x.«»l * - - 

F M = * G « m x4< r o h* 3 - 2 ‘ <?o + 6 x + R x.i> (170) 

L 0n J 


Conversion to vector-dyadic notation and deletion of second-order terms yields 

f x G . " F XJ° + AF x G . > 


where 


and 


% " -G e m x4 /r 0 l P 0 


AF X G = * G e m xA 3 + R X,). 


(171) 


(172) 


(173) 


The resultant gravitational force and torque acting on body X *s obtained by a summation 
over all mass elements contained within the body. 


F«o =‘£f x g o = -G e m x 0- J ? o 


(174) 


i€\ 


.n. 


AF- = CjrfAFO . .G e m x /1- 3 (1 - 30 O 0 O ) • 6, 


Zi 


Jx 00 - **., X fff - 0 


(175) 


(176) 


where 


= ? R x.i * af m = 3 G e^o x *x • X 

*X = ^^ m X4 l»XJ ' R x.«) 1 * R x.i R x.,l 


(177) 


(178) 
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(179) 


For any simulation problem in which gravitational effects must be considered, two 
approaches are available. 

a. The inertially fixed reference is chosen to be fixed at the earth’s center-of-mass. 
Then recall 

Jf, = vector from inertial origin to center of mass of body 1, 


hence 




(180) 


From Equation (154), the vector quantity to be added to row k of |0 E ' j is 


KtS 




X * k, keS L 


k-l 


I?,.,., *?)♦*?! k(S . 


(181) 


2 > 


KtS„ 


k = N + 1 . 


This approach simultaneously solves for both the orbit and the attitude dynamics of the 
modeled spacecraft. It should be recognized, however, that digital solutions may be 
subject to significant numerical error, since the orbit parameters will differ from the 
attitude-dynamics parameters by several orders of magnitude. 

b. The inertially fixed reference is chosen to be fixed at the composite system 
center-of-mass. In this approach, the orbit position vector /f 0 is defined by 
Equations (162) through ( 1 66), and 



(182) 
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Since p 0 is known, F A 0 must be deleted from Equation ( 181) to obtain the 
perturbing force due to gravity-gradient effects alone. Then from Equation ( 1 SI ). 
the vector quantity to be added to row k of |tf> 1 } is 



X = k, keS L 


XeS 


k I 


IW* AF “ +A ^ 


(183) 


keS,, 


k = N + 1. 


THERMALLY INDUCED MOTION 

The effects of appendage deformation due to time-varying thermal gradients can under certain 
conditions adversely influence the attitude dynamics and at times even the attitude stability 
of spacecraft. For example see Reference l ). 

In a directional (solar) thermal field thermal gradients across the diameter of an appendage are 
not established instantaneously, but grow exponentially to a steady state value with a specif- 
ic thermal time constant. The thermal gradient distribution along with the thermal expan- 
sion properties of the appendage define an instantaneous position of thermal equilibrium. 

If the appendage has finite mass and stiffness characteristics the rate at which it will actually 
move to the thermal equilibrium position will be governed by its natural frequencies of 
vibration. 

Virtually all spacecraft attitude dynamics problems attributable to the effects of thermally- 
induced deformation stem from the fact that, relative to the appendage, the direction and 
magnitude of the solar thermal field changes at a system natural frequency. This change in 
thermal input can be caused by such effects as three-axis rotational motion of the space- 
craft. shadowing, or torsional motion of a torsionally weak boom. 

To investigate such problems a crude model of thermal deformation is usually sufficient for 
worst cast type analyses. The particular modeling tool which has been successfully em- 
ployed by the author, within the confines of N-BOD. is a thermal spring; that is. a spring 
which has a time- varying thermal equilibrium position. 


i. 


/ 


S 


52 


Assume that body k has a tendency to thermally deform and that the deformation can 
be adequately modeled as an angular rotation of the body about the free coordinate vector 

Let 


0* M = angular amount body k must rotate at time t about free coordinate 
vector cjf to be in a stale of thermal equilibrium. 

Furthermore let the positic; of thermal equilibrium 0* be governed bv the solution of the 
heat conduction equation 


d 

dt 


e T 



e r 


ST 




(ltw) 


where 


0ST 


> • • •> F ) 


Then if 


thermal time constant; 

a steady state angle of thermal deformation about <f m (calculated 
assuming constant thermal field normal to if ); and 

a function of kinematics which defines the changing magnitude and 
direction of the thermal field relative to the body k fixed reference. 


K m = spring constant of the linear spring which restrains motion about the 
free coordinate vector <f m , 

'^)q m = thermal spring torque about the free coordinate vector cf m which tends 
to drive body k to its position of thermal equilibrium. 

It follows that the net reactive thermal spring torque at ningc point k-1 is 


X. K- 1 ' 2 N (185) 

m @ k-1 

This vector quantity must be added to the k f/ * row of . 
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SYMBOLS 


a Semimajor axis of elliptic orbit 

C x Force associated with the mass of body X and its combined 

centripetal and Coriolis acceleration 

CL m Control torque about free coordinate vector 

D Damping coefficient of linear viscous damper which retards motion 

about free coordinate vector 

*m 



F H 

r k-l 


G 


x 





th] 


Eccentric anomaly 

Eccentricity of elliptic orbit 

Resultant of external forces acting on body k 

Resultant force of constraint acting on body k through hinge point 
k-1 

Resultant force acting on body k due to all mechanisms existing 
between bodies J(k) and k at hinge point k-1 

Resultant external force applied to point i of body k 

Resultant of all external forces applied to body k 

Matrix operator which uses the components of quaternion X £ J(X) to 
define transformation matrix l X '^j( X )L 

Earth’s gravitational constant 

Pseudo-inertia tensor of body X with respect to the origin of nest k - 1 
and the hinge point i - 1 

Linear momentum of body X relative to the inertial origin 
Acceleration of gravity 

Inertial angular momentum of momentum wheel m 

Relative angular momentum of momentum wheel m 

(M X M), square matrix of free coordinate vectors along momentum 
wheel spin axes 


PRECEDING PACK BLANK NOT FILMED 
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Inertia tensor of momentum wheel m about its center-of-mass 

U‘J 


(N + 1 ) X M, rectangular matrix of momentum-wheel inertia tensors 

in 


M X M, square matrix of momentum-wheel inertia tensors 

I(X) 


Body label of the body to which body X is attached at hinge 
point X-l 

K 


Spring constant c ' linear spring which restrains motion about free 
coordinate vector q 

in 

kx 


Inertial angular momentum of body X relative to the inertial origin 

w 

m, m 


Inertial angular momentum of momentum wheel m relative to its 
ccnter-of-mass 



Inertial angular momentum of body X relative to its center-of- 
mass 

M 


Total number of momentum wheels 

MO(m) 


Body label of the gyrostat in which momentum wheel m is embedded 

m x 


Total mass of body X 

m x,i 


Mass of mass element i of body X 

Me 


Mean anomaly 

N 


Total number of rigid bodies, gyrostats, and point masses 

N F 


Total number of free coordinate axes 

n l 


Total number of locked coordinate axes 

N 0 


Total number of coordinate vector triads 

N t 


Total number of free and locked coordinate axes 

r. 


Orbital period 

nx) 


Total number of constrained axes at hinge point X-l of body X 


•k 

Set of unit locked coordinate vectors which span the N L -dimensional 
vector space in which motion is totally constrained 

I?,.?,.- 

• • .?N P ) 

Set of unit free coordinate vectors which span the N F -dimensional 
vector space in which motion is possible 



Quaternion operator, which maps vectors into quaternions 
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LqJ 

Ki 

Sr 

s L 

s 

Ski 

Sr.,, x-i 

&> 

T, 

fx^'c 1 
lx^ ( x)l 
l°x^ ( x)l 


[XI 

(X?,x^x«) 

(X*, Xj, Xj) 


K 


(N + 1 ) X (N F - M), rectangular matrix of free coordinate vectors 

Position vector from center-of-mass of body X to mass point i of 
body X 

Set of all rigid-body, body labels 
Set of all point-mass body labels 
Set of all body labels 

Set of all body labels of those bodies outboard of hinge point k - 1 , 
relative to body 1 . Body labels of those bodies in nest k - 1 

Set of all body labels associated with those bodies lying on the 
topological path from hinge point k - 1 to the center-of-mass of 
body X 

Tensor operator which maps vectors into skew symmetric tensors of 
rank 2, dyads 

Time of perihelion russage 

Transformation matrix, takes vectors from computing frame to 
body X fixed coordinates 

Transformation matrix, takes vectors from body J(X) to body X fixed 
coordinates in the nominal zero stress state 

Transformation matrix, takes vectors from body J(X) to body X fixed 
coordinates at time zero 

True anomaly 

Dyad in *.he k fft row, i th column of [X] 

(N + 1) X (N + 1 ), matrix of inertia and pseudo-inertia dyads 
Coordinate axes defined at inertial origin 

Coordinate axes fixed in body X at hinge point X - 1 

Position vector from hinge point X - 1 of body X to the 
center of mass of body X 

Position vector from hinge point X - 1 to nominal position of 
point mass X 

Position vector from inertial origin to composite N-body system 
center-of-mass 



K 

l 

t 


M 

l»‘l 

li*! 

»„ 

W 

tU 

*» 

* 

k. 


Unit vector aligned with {f 0 

Position vector from inertial origin to hinge point 0 of body 1 

Position vector from hinge point J(X) - 1 to hinge point X - 1 of 
body X 

Position vector from inertial origin to center-of-mass of body X 

Position vector from hinge point k to the center-of-mass of body X 

Skew symmetric tensor form of the vector 7 k x 

Position vector from composite N-body system center-of-mass to the 
center-of-mass of body X 

Rotation quaternion 
Mean motion 

(N + 1 ) X 1 , column matrix of forces and torques associated with 
centripetal and Coriolis acceleration effects 

(N + 1) X 1, column matrix of torques associated with the inertial 
angular momentum of rigid bodies about their own centers-of-mass 

MX 1, column matrix of torques associated with the inertial 
angular momentum of momentum wheels about their own centers- 
of-mass 

Displacement about or along free coordinate vector if 

(N F - M) X 1, column matrix of all 0 m m = 1, 2, . . . , N F - M 

(MX 1), column matrix of all 0 m m = N F - M + 1, . . . , N F 

Inertia tensor of body X about its center-of-mass 

Resultant torque acting on body k due to external causes 

Resultant of external torques acting on body k 

Resultant of external torques acting on momentum wheel m 

Resultant torques of constraint acting on body k through hinge 
point k - 1 



i 


■§;* 


I 


U'*l 

u H 'i 

4M 

i /'I 

-► 

OX 
c w \ 

w x 

M 

)w} 


Resultant torque of constraint acting on momentum wheel m 

Resultant torque acting on body k due to all mechanisms existing 
between bodies J(k) and k at hinge point k - 1 

(N + 1 ) X 1 , column matrix of forces and torques of constraint 
acting between bodies 

M X 1 , column matrix of constraint torques acting on momentum 
wheels 

(N + 1 ) X 1, column matrix of forces and torques due to mechanisms 
acting between bodies 

M X 1 , column matrix of torques due to mechanisms acting on 
momentum wheels 

(N + 1) X 1, column matrix of forces and torques external to total 
system acting on rigid bodies and point masses 

Angular velocity of body X fixed coordinates relative to the 
computing frame fixed coordinates 

Angular velocity of body X fixed coordinates relative to inertially 
fixed coordinates 

Angular velocity of body X fixed coordinates relative to body J(X) 
fixed coordinates 

(N + 1 ) X 1, column matrix of quantities X = 1, . . . , N + 1 
(N + 1 ) X 1 , column matrix of quantities , X = 1 , . . . , N + 1 




Inertial angular velocity of momentum wheel m 

Relative angular velocity of momentum wheel m 

Column matrix 
Square matrix 
Rectangular matrix 
Vector R 

3 X 1 matrix of components of vector R relative to body X fixed 
coordinates 

Quaternion of vector R relative to body X fixed coordinates 
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Time derivative of vector R relative to inertial fixed reference frame 


Time derivative of vector R relative to local reference frame i 

Q 

Same as R when no confusion as to which local reference frame 
differentiation is with respect to 

1 Unit dyad 



£ 

i@k-l 


Summation over all indices i contained in the set S k 


Summation over all indices i of vectors defined at hinge point k - 1 



Multiple product over all indices k contained in set S Q k l in 
decreasing order of magnitude 
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APPENDIX 

QUATERNION TECHNIQUES 


The quaternions of Hamilton constitute a four-dimensional vector space, over the field of 
real numbers, with respect to a basis of four special vectors denoted by 

jl. £,.I 2 .? 3 j . 

The algebraic operations for quaternions are the usual two vector operations of vector 
addition and scalar multiplication, plus the operation of quaternion multiplication. 

In four dimension vector space the quaternion can be defined as 


f * e 0 + e, + e 2 f 2 + e 3 f 3 , 

where 

Ip f 2 > ifjj 

is the set of four linearly independent basis quaternions and 

e o ,e i ,e 2 ,e 3 

are real numbers. The product of any two of the basis quaternions is such that 1 acts as 
the identity and the multiplication table 


8 

= - 1 i= 1,2,3 

(A- 2) 


= "^2^1 ^3 

(A-3) 

Vt 

II 

1 

if 

(A-4) 


* “ ^1 ^3 “ ^2 

(A-5) 


is satisfied. 

Quaternions are most commonly used in conjunction with vectors to describe the effect 
of either a rotation of a vector or a transformation of coordinates. Accordingly, an 
operator is defined which takes vectors into quaternions and vice versa. 


i 
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Let 


and 


If 


£ * quaternion operator which takes vectors into quaternions 


= inverse quaternion operator which maps quaternions into vectors 




J(X) 


J(\) 


column matrix of the components of the vector R relative to 
the body J(X) fixed coordinate axes, 

(A-6) 


then 


*<W JW >- R , W 


= quaternion representation of the vector R given (A-7) 
relative to body J(X) fixed coordinates 


where 

^j(x) = x i ^1 + y 
and 

■*-' <*w ■ 

Euler has shown that an arbitrary rotation of a rigid body about a fixed point is always 
equivalent to a rotation about a line passing through the point. It follows that the relative 
orientation of the coordinate axes fixed in body J(X) and in body X can be completely 
defined by specifying the direction of the line about which and the angle through which 
the coordinate axes { X J , (X \ X* ( *\ X J 3 (X) } must be rotated so that its axes are respectively 
parallel to those of the coordinate axes ] X*, X*, X* } . 

Let 

cf = The unit vector aligned in a right handed sense along the axis about which 
the coordinate axes {x^, X J 2 (X \ X 3 (x) } are rotated, and 


f, * «, f, 


(*!.«• 


(A-8) 


(A-9) 


i 

j 
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I 
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Q = Angle through which the coordinate axes {X , , (x) , X 2 (x) , X 3 (x) | are 
rotated about if so as to be respectively aligned with the coordinate 
axes {x x ,X x ,X x |. 

— ► 

Several texts have shown that the components of R relative to body X fixed coordinates 
can be obtained from the quaternion equation 


where 


and 


R J(X) “ J(X)^X R X J(X)$X 


(A- 10) 

ijx = cos<? / 2 + '^(| < lfj(X)) 

sin 0/2 

(A-l 1) 

JCX}''* = cos ®/2 ~ ^(l^ljcx)) 

| sin 0/2 

(A-l 2) 



(A-13) 


By expressing J(x) f A as 


J(\A " 


+ e 


,1 


+ e, 

2 s 2 


+ e 


3 i 


3 > 


(A- 14) 


it is shown by the direct quaternion multiplication of Equation (A- 10) that 
x i + y Jj + z i h = 

l(ej + e* - e\ - e| ) x + 2 (e, e 2 - e Q e 3 ) y + 2 (c, e 3 + c Q e 2 ) ?.] f, 

+ |2(e 0 e 3 + e, e,) x + (e* - e* + e* - e*) y * 2 (e 2 e 3 - e 0 e,) z| F, 

+ l 2 ( e i e 3 * e o e 2 ) x + 2 ( e 2 e 3 + e o e .) V + < e o ‘ e ? - e l + e 3> Z 1 h, 

(A-l 5) 

f 

i 

i 

I 

v 

? 

k 
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f 



or in matrix notation 




e 2 + e 2 - e 2 - e 2 
e o + e i e 2 e 3 

2 (e. e 2 

- e 0 e 3> 

2 ( e , + e o e 2 

II 

2 ( c o C 3 + e i C 2^ 

e 2 - e 2 
*o e i 

+ e 2 - e 2 

C 2 C 3 

2 ( e 2 e 3 " e 0 e |) 


- 2 ( e ( e 3 - e 0 e 2 ) 

2 ( e 2 e 3 

+ e o e i> 

e l ’ e ? - e l + 


lj(A)' y \l | R }x 


(A- 16) 


where 

= transformation matrix which takes vectors from body X fixed 
coordinates into body J(X) fixed coordinates. 

From this development it can be seen that a matrix operator S can be defined such that 
for the quaternion 


f = e 0 +e,(. + e 2 { 2 + e 3 , 



e 0 + 'l * e l ' e 3 

2(e, e 2 

- e 0 e 3> 

2(e,e 3 

+ e 0 e 2> 

.?(b = 

2( e 0 e 3 + e i tf 2> 

e 2 - e 2 
e o e i 

+ - e 2 

e 2 e 3 

2 ( e 2 e 3 

* e o e i> 


- 2 ( C | e 3 * e 0 e 2> 

2 ( e i e 3 

+ e o e i) 

e 2 - e 2 
e o e i 

■ e 2 + e ; 


Note that from the above definition 


( A- 17) 


(A- 18) 


(A- 19) 


In many practical problems it is infeasible to define the particular eigenvector q and angle 
0 which takes body J(X) fixed axes directly into body X fixed axes. It is possible, however, 
to define several successive rotations through given angles and about given axes which will 
bring the two coordinate axes into alignment. 

Suppose for example that m coordinate rotations are required to conveniently rotate 
body J(X) axes into body X axes. The rotations are as follows: 
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t 


a. Rotate body J(X) axes about cf, thru angle into intermediate reference 
frame 1 , . 

b. Rotate frame I, about cf 2 thru angle 9 J into intermediate reference frame l 2 . 

c. Proceed sequentially until frame I m ., is rotated about q m thru angle 0 m into 
body X fixed frame. 

In quaternion notation this rotation sequence is given by 
R j(\) = J(X)f|, i, f i 2 •• i m . 

where 

, ?, = COS 0 J2 + 

'i-l ‘i 1 

and 

|qj 1 _ components of vector qj relative to the intermediate reference frame 1 M - 

The resultant quaternion which takes body J(X) axes directly into body X axes is 
therefore given by the quaternion product 


** R x i, f * 2 iwff, , (A-20) 

WIT,!,, .,)>««,« 


= J(x)f|, I, f i 2 * (A-22) 

The application of quaternions to the development of coordinate-transformation matrices 
has been demonstrated in the preceeding paragraphs. Quaternions are also used to define 
vector transformation matrices; that is, a quaternion equation can be defined which 
describes the effect of the rotation^ an arbitrary vector R about an eigenvector q 
through the angle 0 into a vector R 1 . If both R and q are defined in the body X fixed 
coordinates, then 


where 


R X = f* \ f 


( A-23) 


f = cos 0/ 2 + J 


(Klx) si 


sin 0/2, 


(A-24) 
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(A-25 1 


^ ( \ ) = |*.)a* 


and j-' ( r a , = {rJ a . <A-2 o) 

■4 

Note that the components of the vector R t arc given relative to the same reference 
frame iri which R and q are given. Furthermore, by making use of the matrix operator 
■f (£) defined by Equations (A- 1 8), Equation ( A-23) can be expressed in matrix notation 
as 






( A-27) 


Direct comparison of Equations (A- 1 6) and ( A-27 ) reveals that the transpose of the 
coordinate transformation matrix is the vector transformation matrix. 
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